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ABSTRACT 

SN 2013dy is a Type la supernova for which we have compiled an extraordinary 
dataset spanning from 0.1 to ~500 days after explosion. We present 10 epochs of 
ultraviolet (UV) through near-infrared (NIR) spectra with HST/STIS, 47 epochs of 
optical spectra (15 of them having high resolution), and more than 500 photometric 
observations in the BVrRiIZYJH bands. SN 2013dy has a broad and slowly declining 
light curve (Ami5(5) = 0.92 mag), shallow SiiiA6355 absorption, and a low velocity 
gradient. We detect strong C II in our earliest spectra, probing unburned progenitor 
material in the outermost layers of the SN ejecta, but this feature fades within a few 
days. The UV continuum of SN 2013dy, which is strongly affected by the metal abun¬ 
dance of the progenitor star, suggests that SN 2013dy had a relatively high-metallicity 
progenitor. Examining one of the largest single set of high-resolution spectra for a 
SN la, we find no evidence of variable absorption from circumstellar material. Com¬ 
bining our UV spectra, NIR photometry, and high-cadence optical photometry, we 
construct a bolometric light curve, showing that SN 2013dy had a maximum lumi¬ 
nosity of lO.Ol^ g X 10^^ erg s“^. We compare the synthetic light curves and spectra 
of several models to SN 2013dy, finding that SN 2013dy is in good agreement with a 
solar-metallicity W7 model. 

Key words: supernovae: general - supernovae: individual (SN 2013dy) 
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1 INTRODUCTION 


et al.|2003 1. A different peak luminosity for SNe la with the 


Type la supernovae (SNe la) are important distance indi¬ 
cators. Observations of SNe la provided the first evidence 
of the accelerating expansion of the universe ( |Riess et al.| 
1998 Perlmutter et al^|1999|. Observational evidence indi¬ 


cates that they are the result of the thermonuclear explo¬ 
sion of an accreting carbon-oxygen white dwarf (WD) star 
in a close binary system (e.g., Hillebrandt fc Niemeyer|2000 


Hillebrandt et aL]|2013 Maoz et al.||2014 |. Although recent 

observations have constrained the exploding star to be con¬ 


sistent with a compact object (Nugent et al. 2011 Bloom 
et al. ][20l2l ), the nature of the presumed companion star is 


not yet clear. 

Some explosion mechanisms have been proposed to ex¬ 
plain the observations of SNe la. For a review of explosion 
models, see Hillebrandt & Niemeyer (20001 and Hillebrandt 


|et al.| ( |2013 1. A pure deflagration model that results from a 
subsonic deflagration flame can produce intermediate-mass 
elements (IMF), but it fails to synthesize enough iron-group 
elements (IGEs) to meet the constraints of normal SNe la. 
However, a pure deflagration may explain some subclasses 


of thermonuclear explosions (Phillips et al. 

2007 Foley et al. 

2013 Kromer et al.||2013 Fink et al.|2014 

. To create more 


IGEs, a model in which a supersonic detonation following 
the deflagration (the so-called delayed-detonation model) is 
proposed. This model produces more IGEs and better ex¬ 
plains normal SNe la than the pure deflagration model (e.g., 
Sim et al.|2013 |. 


Statistical studies using a large sample of SNe la can 
be useful in understanding the physical properties of SN la 
explosions and progenitor systems. However, a well-studied 
single SN la with complete and high-quality observations 
can provide highly constraining information as well. Recent 


SNe la with extensive datasets include SNe 2009ig (Foley 


et al. 

2012a Marion et al.|2013l, 2011fe (e.g., Nugent et al. 

2011 

Li et al.||2011 Bloom et al.||2012l, 2012cg (Silverman 

et al. 

2012b 1, 2012fr ( 

Ghildress et al. 2013b 1 , and 2014J 

(e.g., Zheng et al.|2014| 

Goobar et al. 2014 Amanullah et al. 

2014 Foley et al. 2014 

Marion et al. 1120151. These SNe la 


are all very close {D 20 Mpc), which allows for compre¬ 
hensive datasets including observations at non-optical wave¬ 
lengths and extremely late-time data, and discovered soon 
after explosion, which provides information about the out¬ 
ermost layers of the ejecta. For each of these nearby SNe la, 
large follow-up campaigns were initiated. These high-quality 
datasets have resulted in some of the most important con¬ 
straints for the progenitor systems and explosions of SNe la 
(e.g., progenitor metallicity; Mazzali et al.|[20l4K 


One of the most critical pieces of information one can 
obtain for a SN la is a series of ultraviolet (UV) spectra. 
Theoretical studies show that the metallicity of the progeni¬ 
tor can signihcantly impact both the photometric and spec- 


troscopic properties of SNe la (IHoefiich et al. 

1998 Lentz 

et al. 

|2000| ITimmes et al.||2003| |Kasen et al.||20091 (Walker 

et al. 

20121, but has little effect on optical and near-infrared 


(NIR) spectral properties. While progenitor metallicity does 
not affect the shape of the optical spectral-energy distribu¬ 
tion (SED), the optical spectral features, or the light-curve 
shape dHoefiich et al.|1998 Lentz et al.|200d Podsiadlowski 


et al. 12006 1, it is expected to affect the amount of '^°Ni gener- 

ated in the explosion and thus the peak luminosity (Timmes 


same light-curve shape could introduce a large systematic 
bias in cosmological distance estimates, especially if the av¬ 
erage progenitor metallicity evolves with redshift. 

While observations of SN environments indicate that 


metallicity could subtly affect observables (I Gallagher 

et al. 

2005| (D’Andrea et al.||2011| (Stanishev et al. 2012 Johans- 

son et al.||2013 

Childress et al.|2013a Pan et al.|2014 

1 , the 


most powerful and indicative constraints have come from 


analyses of UV spectra. Foley & Kirshner (20131 examined 


the UV spectra of the “twin” SNe 2011by and 2011fe, which 
had nearly identical optical SEDs, optical colors, and optical 
light-curve shapes. However, their UV spectra differed signif¬ 
icantly, such that the progenitor of SN 2011by had a higher 
metallicity than that of SN 201 Ife. Intriguingly, their peak 
luminosities differed in a way that one would predict based 
on the progenitor metallicity differences; however, some or 
all of this difference may be caused by an incorrect distance 
to SN 2011by ( Graham et al.|2015a |. Regardless of the dis¬ 
tance to SN 2011by, this represents the hrst robust detection 
of different metallicities for SN la progenitors. Additional 
direct analyses of UV spectral time series have placed con¬ 
straints on SN la progenitor metallicity ( [Hachinger et ah] 
2013 Mazzali et al. 12014 1. 

In this work, we present a multi-wavelength analysis of 
SN 2013dy, a SN la known to have the earliest detection 
(~ 2.4 hr after the first light; Zheng et al. 20131. In addi¬ 


tion to high-cadence optical and NIR light curves, over 30 
low-resolution optical spectra, and one of the largest set of 
high-resolution optical spectra of a SN la, our dataset in¬ 
cludes 10 UV spectra observed with the Hubble Space Tele¬ 
scope (HST). We use this high-quality dataset that spans 
epochs from ~0.1 to 500 days after the explosion to con¬ 
strain the properties of the explosion and progenitor system 
for SN 2013dy. 

This manuscript is organized as follows. In § [^ we 
present the photometric and spectroscopic observations of 
SN 2013dy. In § [^ and we present the photometric 
and spectroscopic properties of SN 2013dy, respectively. 
We compare models to our data in § |^ and summarize 
our hndings in § [^ Throughout this paper, we assume 
Ho = 70kms“^ Mpc“^ and a flat universe with Hm = 0.3. 


2 OBSERVATIONS 

SN 2013dy was discovered in NGC 7250 {D — 13.7 Mpc; 


Tully et al. 20091 on 2013 July 10.46 (UT dates are used 


throughout this paper) and classihed as a young SN la 
( Gasper et al.|[2M3 l. Radio observations taken ~1 week af¬ 
ter maximum brightness resulted in nondetections (Perez- 
Torres et al.|2013|. 


Zheng et al. (20131 obtained early-time photometry 
of SN 2013dy starting ~0.1d after explosion. These data 
constrained the progenitor star to be a compact object 
{Ro ^ 0.25 Rq). Their earliest spectrum had very strong 
C II lines, tracing the unburned material in the outer ejecta. 
We combine the photometric and spectroscopic data from 


Zheng et al. (20131 with our own data presented below. All 


of the data will be made available in the WISeREP archive 
( Yaron fc Gal-Yam|2012 I. 
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Figure 1. The KAIT BVRI and RATIR riZYJH light curves 
of SN 2013dy. The solid lines represent the best fit (from SNooPy) 
to the observed light curves. 


2.1 Photometry 


We obtained broadband BVRI photometry of SN 2013dy 
with the 0.76-m Katzman Automatic Imaging Telescope 
(KAIT; Filippenko et al. 2001| . The multi-band images were 
observed with the KAIT4 filter set from —16 d to -1-337 d rel¬ 
ative to B-band maximum (MJD = 56501.105). The data 


were reduced using our image-reduction pipeline (Gane- 


shalingam et al. 20101. The point-spread function (PSF) 
photometry was performed using DOAPHOT (Stetson 19871. 


The SN instrumental magnitudes were calibrated to nearby 
Sloan Digital Sky Survey (SDSS) standard stars and trans¬ 
formed into the Landolt system. Details about the KAIT 
photometry are provided by Zheng et al. (20131. In Table 
and we list the KAIT photometry used in this work. 

We also obtained riZYJH photometry of SN 2013dy 
with the multi-channel Reionization And Transients In¬ 


fraRed camera (RATIR; Butler et al. 20121 mounted on 
the 1.5-m Johnson telescope at the Mexican Observatorio 
Astronorhico Nacional on Sierra San Pedro Martir in Baja 


California, Mexico (Watson et al. 20121. Typical observa¬ 


tions include a series of 80-s exposures in the ri bands and 
60-s exposures in the ZY JH bands, with dithering between 
exposures. RATIR’s fixed IR filters cover half of their re¬ 
spective detectors, automatically providing off-target IR sky 
exposures while the target is observed in the neighboring fil¬ 


ter. Master IR sky frames are created from a median stack of 
off-target images in each IR filter. No off-target sky frames 
were obtained on the optical CCDs, but the small galaxy 
size and sufficient dithering allowed for a sky frame to be 
created from a median stack of all the images in each filter. 
Flat-field frames consist of evening sky exposures. Given the 
lack of a cold shutter in RATIR’s design, IR dark frames are 
not available. Laboratory testing, however, confirms that the 


dark current is negligible in both IR detectors (Fox et al. 


20121 . 


The RATIR data were reduced, coadded, and anal¬ 
ysed using standard CCD and IR processing techniques in 
IDL and Python, utilizing the online astrom etry programs 
SExtractor and SWarj0(s 


Foley et al. 


2014 


and references 

therein). Calibration was performed using field stars with re¬ 
ported fluxes in both 2MASS ( |Skrutskie et al.|2006| and the 
SDSS Data Release 9 Catalogue ( Ahn et al.||2012 |. Table 
lists the RATIR photometry. In Figure [fj we present both 
the KAIT and RATIR light curves and the best-fit template 
light curves (see § 3.11. 


2.2 Spectroscopy 


As part of our follow-up campaign, we obtained a series 
of low-resolution near-UV through NIR spectra with HST, 
as well as low- and high-resolution optical spectra from a 
variety of sources. 

We procured 10 epochs of near-UV through NIR spec¬ 
troscopy with the Space Telescope Imaging Spectrograph 
(STIS) on HST. For each epoch, we obtained data with three 
different setups to cover the entire wavelength range: the 
NUV MAMA detector with the G230L grating, the CCD 
detector with the G430L grating, and the CCD detector 
with the G750L grating. Combined, the three setups have 
a wavelength range of 1600-10,230 A. The HST data were 
reduced using the standard Space Telescope Science Data 
Analysis System (STSDAS) routines, including bias subtrac¬ 
tion, flat-fielding, wavelength calibration, and flux calibra¬ 
tion ( Foley et al.|2012b (. The HST observations correspond 
to phases ranging from —7 to -1-21 d relative to B-band max¬ 
imum brightness. A log of our HST spectroscopic observa¬ 
tions is presented in Table and the spectra are shown in 
Figure 

In addition to the HST spectroscopy, we obtained 32 
optical spectra of SN 2013dy with phases of —16 to -1-480 d 
relative to B-band maximum brightness with a variety of 
ground-based facilities. The 8 earliest spectra in the dataset 


were already studied by Zheng et al. (20131, but are included 


in our analysis for completeness. The new spectra (24 of the 
32 total) were observed with the Low Resolution Imaging 
Spectrometer (LRIS; Oke et al.|1995 1 mounted on the 10-m 
Keck-I telescope, the DEep Imaging Multi-Object Spectro¬ 
graph (DEIMOS; [Faber et al. 20031 on the 10-m Keck-II 
telescope, the Kast Spectrograph ( Miller fc Stone||l993 ( on 
the Lick Observatory 3-m Shane telescope, and the Low- 


Resolution Spectrograph (LRS; Hill et al^|1998 1 on the 9.2-m 
Hobby-Eberly Telescope (HET). 

Standard CCD processing and spectrum extraction 


^ SExtractor and SWarp can be accessed from 
http: / / www.astromatic.net / software 
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Table 1. KAIT photometry of SN 2013dy. 


MJD 

B 

V 

R 

I 


(mag) 

(mag) 

(mag) 

(mag) 


56484.42 

16.76 

± 

0.04 

16.41 

± 

0.03 

16.17 

± 

0.03 

16.07 

± 

0.04 

56485.46 

16.21 

± 

0.03 

15.84 

± 

0.03 

15.60 

± 

0.03 

15.46 

± 

0.03 

56486.45 

15.77 

± 

0.07 

15.37 

± 

0.05 

15.17 

± 

0.03 

15.03 

± 

0.03 

56487.38 

15.33 

± 

0.03 

14.96 

± 

0.02 

14.74 

± 

0.02 

14.62 

± 

0.02 

56488.40 

14.88 

± 

0.03 

14.56 

± 

0.03 

14.34 

± 

0.03 

14.23 

± 

0.04 

56489.37 

14.55 

± 

0.03 

14.24 

± 

0.03 

14.02 

± 

0.03 

13.91 

± 

0.04 

56490.43 

14.24 

± 

0.02 

13.95 

± 

0.02 

13.72 

± 

0.02 

13.63 

± 

0.02 

56491.40 

14.04 

± 

0.03 

13.74 

± 

0.02 

13.53 

± 

0.03 

13.44 

± 

0.04 

56492.39 

13.86 

± 

0.03 

13.57 

± 

0.03 

13.35 

± 

0.03 

13.29 

± 

0.03 

56493.40 

13.73 

± 

0.06 

13.42 

± 

0.04 

13.22 

± 

0.03 

13.17 

± 

0.04 

56494.41 

13.61 

± 

0.03 

13.31 

± 

0.03 

13.11 

± 

0.03 

13.07 

± 

0.03 

56496.32 

13.46 

± 

0.05 

13.15 

± 

0.03 

12.98 

± 

0.02 

12.97 

± 

0.02 

56497.41 

13.37 

± 

0.03 

13.06 

± 

0.02 

12.93 

± 

0.02 

12.95 

± 

0.03 

56498.38 

13.33 

± 

0.03 

13.03 

± 

0.02 

12.90 

± 

0.02 

12.96 

± 

0.02 

56499.38 

13.33 

± 

0.07 

12.99 

± 

0.05 

12.87 

± 

0.06 

12.96 

± 

0.06 

56500.45 

13.28 

± 

0.03 

12.96 

± 

0.03 

12.83 

± 

0.03 

12.97 

± 

0.03 

56501.39 

13.30 

± 

0.03 

12.96 

± 

0.03 

12.84 

± 

0.03 

13.00 

± 

0.04 

56502.35 

13.30 

± 

0.03 

12.96 

± 

0.03 

12.83 

± 

0.03 

13.02 

± 

0.03 

56503.30 

13.33 

± 

0.03 

12.96 

± 

0.03 

12.84 

± 

0.03 

13.06 

± 

0.03 

56504.35 

13.37 

± 

0.03 

12.98 

± 

0.03 

12.84 

± 

0.03 

13.09 

± 

0.03 

56505.34 

13.40 

± 

0.03 

13.00 

± 

0.02 

12.88 

± 

0.02 

13.14 

± 

0.02 

56506.35 

13.44 

± 

0.02 

13.03 

± 

0.02 

12.90 

± 

0.02 

13.19 

± 

0.02 

56507.33 

13.50 

± 

0.04 

13.07 

± 

0.03 

12.93 

± 

0.04 

13.24 

± 

0.04 

56508.33 

13.54 

± 

0.02 

13.09 

± 

0.02 

12.99 

± 

0.02 

13.30 

± 

0.03 

56509.34 

13.59 

± 

0.03 

13.13 

± 

0.03 

13.05 

± 

0.03 

13.38 

± 

0.03 

56510.34 

13.66 

± 

0.03 

13.18 

± 

0.02 

13.11 

± 

0.03 

13.43 

± 

0.03 

56511.31 

13.73 

± 

0.03 

13.21 

± 

0.02 

13.15 

± 

0.02 

13.48 

± 

0.02 

56512.31 

13.82 

± 

0.03 

13.30 

± 

0.02 

13.26 

± 

0.03 

13.59 

± 

0.03 

56513.32 

13.90 

± 

0.03 

13.34 

± 

0.03 

13.33 

± 

0.02 

13.64 

± 

0.03 

56514.32 

14.01 

± 

0.03 

13.42 

± 

0.03 

13.39 

± 

0.03 

13.69 

± 

0.04 

56516.33 

14.19 

± 

0.03 

13.54 

± 

0.03 

13.50 

± 

0.03 

13.73 

± 

0.03 

56517.32 

14.30 

± 

0.04 

13.60 

± 

0.03 

13.54 

± 

0.03 

13.72 

± 

0.03 

56518.36 

14.39 

± 

0.03 

13.65 

± 

0.02 

13.56 

± 

0.02 

13.72 

± 

0.02 

56519.31 

14.51 

± 

0.03 

13.72 

± 

0.03 

13.60 

± 

0.02 

13.72 

± 

0.02 

56520.34 

14.60 

± 

0.03 

13.75 

± 

0.03 

13.59 

± 

0.02 

13.66 

± 

0.03 

56521.41 

14.74 

± 

0.15 

13.80 

± 

0.11 

13.60 

± 

0.06 

13.62 

± 

0.08 

56522.33 

14.80 

± 

0.03 

13.86 

± 

0.02 

13.63 

± 

0.02 

13.64 

± 

0.02 

56523.25 

14.90 

± 

0.03 

13.89 

± 

0.02 

13.64 

± 

0.02 

13.62 

± 

0.02 

56524.31 

14.96 

± 

0.07 

13.95 

± 

0.05 

13.64 

± 

0.05 

13.57 

± 

0.04 

56525.29 

15.07 

± 

0.04 

13.99 

± 

0.03 

13.68 

± 

0.03 

13.58 

± 

0.03 

56527.35 

15.24 

± 

0.04 

14.09 

± 

0.03 

13.71 

± 

0.03 

13.56 

± 

0.03 

56528.36 

15.32 

± 

0.03 

14.12 

± 

0.02 

13.72 

± 

0.03 

13.54 

± 

0.03 

56529.33 

15.40 

± 

0.03 

14.15 

± 

0.03 

13.74 

± 

0.03 

13.51 

± 

0.03 

56530.33 

15.48 

± 

0.03 

14.22 

± 

0.02 

13.78 

± 

0.02 

13.52 

± 

0.03 

56531.32 

15.54 

± 

0.03 

14.25 

± 

0.03 

13.79 

± 

0.03 

13.50 

± 

0.03 

56532.35 

15.62 

± 

0.03 

14.31 

± 

0.02 

13.84 

± 

0.02 

13.51 

± 

0.02 

56533.30 

15.67 

± 

0.03 

14.35 

± 

0.03 

13.87 

± 

0.03 

13.51 

± 

0.03 

56534.31 

15.73 

± 

0.03 

14.41 

± 

0.02 

13.93 

± 

0.02 

13.55 

± 

0.02 

56535.28 

15.80 

± 

0.03 

14.46 

± 

0.02 

13.98 

± 

0.02 

13.57 

± 

0.03 

56536.26 

15.85 

± 

0.03 

14.50 

± 

0.02 

14.02 

± 

0.02 

13.58 

± 

0.03 

56537.29 




14.63 

± 

0.10 







56538.28 




14.76 

± 

0.12 

14.09 

± 

0.15 

13.63 

± 

0.09 

56539.28 

16.01 

± 

0.03 

14.69 

± 

0.03 

14.21 

± 

0.03 

13.79 

± 

0.03 

56540.29 

16.00 

± 

0.03 

14.70 

± 

0.02 

14.25 

± 

0.02 

13.85 

± 

0.02 

56541.26 

16.04 

± 

0.04 

14.77 

± 

0.03 

14.34 

± 

0.03 

13.93 

± 

0.03 

56542.28 

16.10 

± 

0.03 

14.83 

± 

0.02 

14.38 

± 

0.02 

13.99 

± 

0.02 

56543.26 

16.13 

± 

0.03 

14.86 

± 

0.02 

14.43 

± 

0.02 

14.04 

± 

0.03 

56544.24 

16.17 

± 

0.03 

14.90 

± 

0.03 

14.47 

± 

0.03 

14.10 

± 

0.03 

56545.26 

16.20 

± 

0.04 

14.94 

± 

0.02 

14.54 

± 

0.02 

14.18 

± 

0.02 

56546.26 

16.19 

± 

0.04 

14.97 

± 

0.02 

14.56 

± 

0.03 

14.22 

± 

0.03 

56548.24 

16.24 

± 

0.03 

15.03 

± 

0.02 

14.64 

± 

0.02 

14.33 

± 

0.02 

56551.26 

16.28 

± 

0.04 

15.13 

± 

0.03 

14.75 

± 

0.03 

14.48 

± 

0.03 

56553.24 

16.38 

± 

0.06 

15.21 

± 

0.03 

14.84 

± 

0.02 

14.60 

± 

0.03 

56555.27 

16.34 

± 

0.05 

15.23 

± 

0.03 

14.88 

± 

0.03 

14.67 

± 

0.03 

56558.23 

16.41 

± 

0.04 

15.32 

± 

0.02 

14.98 

± 

0.02 

14.82 

± 

0.03 

56561.23 

16.49 

± 

0.04 

15.41 

± 

0.02 

15.09 

± 

0.02 

14.98 

± 

0.02 

56563.22 

16.45 

± 

0.03 

15.44 

± 

0.02 

15.14 

± 

0.03 

15.04 

± 

0.03 

56565.22 

16.49 

± 

0.03 

15.50 

± 

0.03 

15.21 

± 

0.03 

15.13 

± 

0.03 

56567.23 

16.54 

± 

0.06 

15.60 

± 

0.03 

15.28 

± 

0.03 

15.23 

± 

0.04 
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Table 2. KAIT photometry of SN 2013dy (continued). 


MJD 

B 

V 

R 

/ 


(mag) 

(mag) 

(mag) 

(mag) 


56569.19 

16.56 

± 

0.05 

15.61 

± 

0.03 

15.37 

± 

0.02 

15.35 

± 

0.03 

56571.20 

16.58 

± 

0.04 

15.67 

± 

0.03 

15.39 

± 

0.03 

15.40 

± 

0.03 

56573.20 

16.59 

± 

0.05 

15.70 

± 

0.03 

15.45 

± 

0.03 

15.49 

± 

0.04 

56575.21 

16.68 

± 

0.05 

15.79 

± 

0.02 

15.52 

± 

0.02 

15.60 

± 

0.03 

56577.20 

16.64 

± 

0.04 

15.81 

± 

0.03 

15.56 

± 

0.03 

15.66 

± 

0.03 

56579.18 

16.69 

± 

0.05 

15.89 

± 

0.03 

15.63 

± 

0.03 

15.73 

± 

0.03 

56581.18 

16.73 

± 

0.06 

15.93 

± 

0.03 

15.71 

± 

0.03 

15.84 

± 

0.04 

56583.19 

16.79 

± 

0.07 

15.96 

± 

0.03 

15.75 

± 

0.03 

15.89 

± 

0.04 

56585.19 

16.81 

± 

0.06 

16.00 

± 

0.03 

15.82 

± 

0.03 

15.98 

± 

0.03 

56587.18 

16.82 

± 

0.04 

16.08 

± 

0.04 

15.88 

± 

0.03 

16.05 

± 

0.04 

56589.18 

16.88 

± 

0.04 

16.11 

± 

0.03 

15.93 

± 

0.03 

16.13 

± 

0.04 

56592.16 

16.91 

± 

0.05 

16.19 

± 

0.04 

16.01 

± 

0.03 

16.24 

± 

0.03 

56596.19 

17.07 

± 

0.05 

16.28 

± 

0.03 

16.15 

± 

0.03 

16.41 

± 

0.04 

56598.14 

17.00 

± 

0.04 

16.32 

± 

0.03 

16.19 

± 

0.03 

16.48 

± 

0.04 

56600.14 

17.02 

± 

0.16 

16.40 

± 

0.12 

16.24 

± 

0.21 

16.88 

± 

0.14 

56602.13 

17.04 

± 

0.04 

16.42 

± 

0.03 

16.32 

± 

0.03 

16.57 

± 

0.04 

56604.13 

17.08 

± 

0.05 

16.49 

± 

0.03 

16.38 

± 

0.03 

16.69 

± 

0.05 

56606.13 

17.14 

± 

0.05 

16.52 

± 

0.04 

16.43 

± 

0.04 

16.73 

± 

0.05 

56609.12 

17.14 

± 

0.07 

16.56 

± 

0.04 

16.49 

± 

0.05 

16.73 

± 

0.06 

56614.14 

17.19 

± 

0.09 

16.65 

± 

0.06 

16.61 

± 

0.05 

16.88 

± 

0.10 

56618.09 

17.29 

± 

0.08 

16.78 

± 

0.04 

16.78 

± 

0.04 

17.11 

± 

0.08 

56639.13 

17.67 

± 

0.02 

17.15 

± 

0.02 

17.33 

± 

0.02 

17.47 

± 

0.04 

56642.13 

17.73 

± 

0.06 

17.24 

± 

0.03 

17.53 

± 

0.08 

17.55 

± 

0.09 

56644.10 




17.39 

± 

0.12 

17.46 

± 

0.13 




56646.07 




17.35 

± 

0.10 

17.45 

± 

0.09 




56647.10 

17.84 

± 

0.14 

17.38 

± 

0.08 

17.65 

± 

0.11 




56656.11 

17.92 

± 

0.10 

17.58 

± 

0.07 

17.75 

± 

0.07 

18.00 

± 

0.13 

56665.11 

18.13 

± 

0.13 

17.73 

± 

0.08 

17.92 

± 

0.13 

18.20 

± 

0.18 

56665.12 

18.10 

± 

0.03 

17.65 

± 

0.04 

18.02 

± 

0.05 

17.97 

± 

0.09 

56667.15 

18.13 

± 

0.10 

17.69 

± 

0.06 







56675.12 

18.20 

± 

0.02 

17.77 

± 

0.02 

18.15 

± 

0.03 

18.12 

± 

0.05 

56676.10 

18.29 

± 

0.25 

17.85 

± 

0.15 

18.04 

± 

0.18 

18.43 

± 

0.37 

56678.11 

18.26 

± 

0.02 

17.82 

± 

0.02 

18.24 

± 

0.03 

18.13 

± 

0.06 

56789.46 

20.21 

± 

0.28 










56838.39 

20.84 

± 

0.24 

20.22 

± 

0.19 

21.13 

± 

0.59 

20.03 

± 

0.28 


Table 3. RATIR photometry of SN 2013dy. 


MJD 

r 

(mag) 

i 

(mag) 

Z 

(mag) 

Y 

(mag) 

J 

(mag) 

H 

(mag) 

56491.30 

13.86 ±0.02 

13.96 ± 

0.02 

13.93 ±0.02 

14.15 ± 

0.02 

14.43 ±0.05 

14.93 ± 0.07 

56498.28 

13.01 ±0.02 

13.32 ± 

0.02 

13.26 ±0.02 

13.66 ± 

0.02 

13.79 ±0.05 

14.39 ± 0.07 

56501.27 

12.89 ±0.02 

13.35 ± 

0.02 

13.22 ±0.02 

13.89 ± 

0.02 



56503.30 

12.90 ±0.02 

13.40 ± 

0.02 

13.18 ±0.02 

14.10 ± 

0.02 

14.02 ±0.05 

14.64 ±0.07 

56504.33 

12.88 ±0.02 

13.44 ± 

0.02 

13.20 ±0.02 

14.20 ± 

0.02 

14.09 ±0.05 

14.63 ± 0.07 

56505.45 

12.93 ±0.02 

13.48 ± 

0.02 

13.24 ±0.02 

14.29 ± 

0.02 

14.19 ±0.05 

14.69 ± 0.07 

56506.45 

12.95 ±0.02 

13.52 ± 

0.02 

13.29 ±0.02 

14.35 ± 

0.02 

14.28 ±0.05 

14.68 ± 0.07 

56508.32 

13.03 ±0.02 

13.62 ± 

0.02 

13.34 ±0.02 

14.45 ± 

0.02 

14.55 ±0.05 

14.73 ± 0.07 

56509.44 

13.09 ±0.02 

13.69 ± 

0.02 

13.40 ±0.02 

14.51 ± 

0.02 



56510.43 

13.12 ±0.02 

13.76 ± 

0.02 

13.44 ±0.02 

14.54 ± 

0.02 

14.89 ±0.05 

14.76 ± 0.07 

56511.35 

13.21 ±0.02 

13.83 ± 

0.02 

13.46 ±0.02 

14.57 ± 

0.02 

15.05 ±0.05 

14.79 ± 0.07 

56512.39 

13.30 ±0.02 

13.93 ± 

0.02 

13.49 ±0.02 

14.57 ± 

0.02 

15.26 ±0.05 

14.81 ± 0.05 

56516.26 

13.53 ±0.02 

14.19 ± 

0.02 

13.55 ±0.02 

14.49 ± 

0.02 

15.66 ±0.05 

14.74 ±0.05 

56517.24 

13.56 ±0.02 

14.19 ± 

0.02 

13.54 ±0.02 

14.44 ± 

0.02 

15.64 ±0.05 

14.72 ±0.05 

56519.47 

13.63 ±0.02 

14.15 ± 

0.02 

13.49 ±0.02 

14.31 ± 

0.02 

15.62 ±0.05 

14.63 ± 0.05 

56521.42 

13.66 ±0.02 

14.10 ± 

0.02 

13.47 ±0.02 

14.17 ± 

0.02 

15.55 ±0.05 

14.57 ±0.05 

56522.30 

13.67 ±0.02 

14.09 ± 

0.02 

13.47 ±0.02 

14.12 ± 

0.02 

15.53 ±0.05 

14.55 ±0.05 

56535.30 


14.04 ± 

0.02 

13.46 ±0.05 





56541.30 

14.39 ±0.02 

14.36 ± 

0.02 

14.00 ±0.02 

13.62 ± 

0.02 

15.39 ±0.05 

14.75 ±0.05 

56545.30 









56552.30 

14.82 ±0.02 

14.89 ± 

0.02 

14.68 ±0.02 

14.25 ± 

0.02 

16.17 ±0.05 

15.35 ±0.05 
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Figure 2. Left: Spectra of SN 2013dy taken with HST/STIS, with arbitrary vertical offsets for clarity. Right: Same as left panel, but a 
closer look at the UV region. The spectra smoothed with a Savitzky-Golay filter ||Savitzky fc Golay|1964ll are shown in red curves. 


Table 4. Log of HST spectroscopic observations of SN 2013dy. 


Date 

(MJD) 

Phase 

(days) 

Exposure time 

(s) 




G230L 

G420L G750L 

56494.48 

-6.6 

3618 

336 

336 

56498.60 

-2.5 

1382 

64 

64 

56500.32 

-0.8 

1382 

64 

64 

56502.31 

-fl.2 

1382 

64 

64 

56505.57 

-1-4.5 

1382 

64 

64 

56509.49 

-1-8.4 

1382 

64 

64 

56513.09 

-1-12.0 

1382 

64 

64 

56515.14 

-1-14.0 

1382 

64 

64 

56518.99 

-1-17.9 

1382 

64 

64 

56521.91 

-1-20.8 

3618 

336 

336 


were accomplished with IRAF. The data were extracted us¬ 
ing the optimal algorithm of Horne (19861. Low-order poly¬ 
nomial fits to calibration-lamp spectra were used to establish 
the wavelength scale, and small adjustments derived from 
night-sky lines in the object frames were applied. We em¬ 
ployed our own IDL routines to flux calibrate the data and 
remove telluric lines using the well-exposed continua of the 


spectrophotometric standard stars (Wade & Horne 1988 Fo- 
|ley et al.|2003[ ). Details of our spectroscopic reduction tech¬ 
niques are described by Silverman et al. ( 2012a[ ). A log of 
spectroscopic observations is presented in Table The com¬ 
plete spectral sequence can be found in Figure!^ 

We also observed SN 2013dy with various high- 
resolution spectrographs on 15 different occasions spanning 


Table 5. Log of low-resolution optical spectroscopic observations 
of SN 2013dy. 


Date 

(MJD) 

Phase 

(day) 

Instrument 

Exp. time 

(s) 

56497.28 

-3.8 

HET/LRS 

380 

56503.26 

-1-2.2 

HET/LRS 

577 

56506.24 

-b5.1 

HET/LRS 

100 

56506.44 

-1-5.3 

Keck/DEIMOS 

60 

56508.22 

-f7.1 

HET/LRS 

260 

56508.41 

-1-7.3 

Lick/Kast 

360 

56511.21 

-1-10.1 

HET/LRS 

400 

56512.41 

-1-11.3 

Lick/Kast 

360 

56515.45 

-1-14.3 

HET/LRS 

200 

56516.20 

-1-15.1 

HET/LRS 

200 

56516.41 

-1-15.3 

Lick/Kast 

360 

56534.42 

-1-33.3 

Lick/Kast 

180 

56541.32 

-1-40.2 

Lick/Kast 

180 

56545.36 

-1-44.3 

Keck/DEIMOS 

300 

56570.23 

-1-69.1 

Lick/Kast 

450 

56575.26 

-1-74.1 

Lick/Kast 

450 

56591.33 

-1-90.2 

Lick/Kast 

900 

56599.24 

-1-98.1 

Lick/Kast 

1800 

56624.24 

-1-123.1 

Lick/Kast 

1800 

56629.22 

-1-127.9 

Keck/LRIS 

350 

56632.12 

+131.0 

Lick/Kast 

1800 

56834.58 

+333.5 

Keck/DEIMOS 

2400 

56924.34 

+423.0 

Keck/LRIS 

1000 

56981.32 

+480.0 

Keck/LRIS 

1000 
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Figure 3. The optical spectra of SN 2013dy. 


Table 6. Log of high-resolution observations of SN 2013dy. 


Date 

(MJD) 

Phase 

(day) 

Instrument 

Texp. 

(s) 

Resolution 

56491.42 

-9.7 

T illinghast/TRES 

2831 

30,000 

56494.09 

-7.0 

Mercator/HERMES 

3600 

85,000 

56494.26 

-6.8 

HET/HRS 

1000 

33,600 

56495.13 

-6.0 

Mercator/HERMES 

3600 

85,000 

56497.26 

-3.8 

HET/HRS 

1000 

33,600 

56502.23 

-fl.l 

VLT/UVES 

600 

40,000 

56504.25 

+3.1 

HET/HRS 

1000 

33,600 

56505.05 

-1-3.9 

Mercator/HERMES 

2700 

85,000 

56506.61 

-1-5.5 

Keck/HIRES 

300 

48,000 

56511.06 

-1-10.0 

Mercator/HERMES 

2700 

85,000 

56513.22 

-1-12.1 

HET/HRS 

1200 

33,600 

56515.97 

-1-14.9 

TNG/HARPS 

1800 

115,000 

56516.20 

-1-15.1 

HET/HRS 

1200 

33,600 

56519.08 

-1-18.0 

Mercator/HERMES 

3600 

85,000 

56520.07 

-1-19.0 

Mercator/HERMES 

3600 

85,000 


for the remaining spectra, we produced synthetic telluric 
spectra ( [Smette et al.|[2015[ ). We used B-splines to fit the 
continuum and normalize the spectra. 


3 PHOTOMETRIC ANALYSIS 

In this section, we analyze the SN 2013dy photometry. We fit 
the light curves, deriving estimates of the host-galaxy red¬ 
dening and extinction, examine the color evolution, and con¬ 
struct a bolometric light curve. Using broadband photome¬ 
try of the host galaxy, we derive host parameters including 
stellar mass (Msteiiar) and star-formation rate (SFR). 


a phase range of —10 to -1-19 days. This corresponds to a 
mean (and median) cadence of 2 days, with the largest gap 
between spectra being 5 days. The high-resolution spectra 
were obtained with the Tillinghast Reflector Echelle Spec¬ 
trograph (TRES) mounted on the 1.5-m Tillinghast tele¬ 
scope, the High Efficiency and Resolution Mercator Echelle 
Spectrograph (HERMES; [Raskin et ah 20111 on the 1.2- 
m Mercator Telescope, the High-Resolution Spectrograph 


(HRS; Tull 19981 on the HET, the Ultraviolet and Visual 
Echelle Spectrograph (UVES; Dekker et al. 20001 on the 
Very Large Telescope (VLT), the High Resolution Echelle 
Spectrometer (HRES; Vogt et al. 19941 on the Keck-I 
telescope, and the High Accuracy Radial velocity Planet 
Searcher (HARPS-N; Cosentino et al.|2012 1 on the Telescope 
Nationale Galileo (TNG). These data represent one of the 
largest single sets of high-resolution spectra of a SN The 
spectra have resolutions ranging from ~ 30,000 to 115,000. A 
log of observations is given in Table The spectra were re- 
duced with standard procedures (e.g., [Sternberg et al.|2014 


and references therein) including telluric correction. For a 
subset of the observations, we observed standard stars to 
determine the telluric absorption in the SN 2013dy spectra; 


^ There are several separate sets of high-resolution spectra of 
SN 2 014J ||Foley et al.||2014| [Goobar et ^ |2014[ [Welty et al.| 


|2014| [Graham et al.||20i5b| | Jack et al.||2015 l; the largest single 
set of spectra has 33 epochs. The largest single set of spectra of 
SN 201 Ife has 12 epochs I Patat et al.|20i3l. 


3.1 Light-Curve Fitting 

We fit the SN 2013dy BVrRiIZYJH light curves simultane¬ 
ously with SNooPy (Burns et al.[[2011|. SNooPy is a Python- 
based light-curve fitter, extending the method of |Prieto[ 
et al. (20061. The decline-rate parameter. Amis, (similar 
to Ami 5 (B) defined by Phillips] 1993 corresponding to the 
R-band decline 15 days after maximum brightness) is used 
to parameterize the SN light-curve shape in SNooPy. Given 
that our dataset were observed with different photometric 
systems, we provide SNooPy with corresponding filter func¬ 
tions and photometric zero-points to perform K- and S- 
corrections. We adopted the default “EBV model” in SNooPy 
to fit the light curves, which is described in detail below. 
The fitting results in measurements of Amis, time of B- 
band maximum brightness, distance modulus (DM) and the 
host-galaxy reddening E{B — U)host. 

For the EBV model, a R-band peak absolute magni¬ 
tude (Mb) as well as the colors are assumed based on the 
value of Amis. Six different calibrations are provided in the 
model. The parameters of these calibrations were derived by 
Folatelli et al.| ( 2010| t o minimize the Phillips relation (Mb 
versus Amis; Phillips 19931 using different subsets of the 


training SN sample, with different Rv determined (or fixed) 
for each calibration. 

We tested all six available calibrations and found that 
calibration #2 (see Table 9 in Folatelli et al.|2(il0 1 with the 
EBV model, which corresponds to Rv = 3.1, produces the 
best fitting of the SN 2013dy light curves (in terms of 
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Table 7. Results of SNooPy light-curve fitting for SN 2013dy“. 


Tmax 

(MJD) 

-^max 

(mag) 

Amis 

(mag) 

DM 

(mag) 

E{B - y)host 

(mag) 

56501.105 

13.229 ±0.010 

0.886 ± 0.006 

31.488 ±0.010 

0.206 ± 0.005 


“The statistical uncertainties of the measurements are generated by Monte-Carlo simulations. 


with Xi/ = 2.93). We adopted this setup when fitting the 
light curves. The results of the light-curve fitting are shown 
in Table The light-curve fitting with different reddening 
models gives consistent Amis, but lower E{B — V)host for 
those with low Rv {E{B — F)host = 0.16 and 0.18 mag for 
Rv = 1.46 and 1.01, respectively). This leads to a lower 
extinction measurement compared to that derived from a 
Rv = 3.1 reddening model. However, the low-i?v models 
produce much poorer fits to the light curves, especially in 
the NIR, than Rv = 3.1. The setnp with Rv = 3.1 is 
still strongly preferred for SN 2013dy even if we only fit 
the BVRI light curves (where the data were better sam¬ 
pled). Previous studies have found that a low Rv value is 
preferred for some highly-reddened SNe la (e.g., |Krisciunas| 
et al.|[2006 Wang et al.||2008 l. On the other hand, we find 


that the dust reddening of SN 2013dy is much more likely 
to have a “high” value of Rv {Rv ~ 3.1) than it is a “low” 
value {Rv < 2). For the rest of the study, we therefore use 
the EBV-i?v =3.1 model for our primary results, but ex¬ 
amine how the conclusions change with a lower value of Rv 
when appropriate. 

SNooPy measures a decline-rate parameter Amis = 
0.89. This template-derived Amis is similar to the conven¬ 
tional B-band decline rate Amis(B), but with some sys¬ 
tematic differences. Using the relation provided by |Burns| 
[er5:] ( (MT| , we calculate a B-band decline rate Ami 5 (B) = 
0.92 mag. This indicates SN 2013dy has a relatively slower 
decline rate, compared to the normal SN la like SN 2011fe 
(Amis(B) « 1.1 mag; Munari et al. 2013| ), but is similar 
to another slow decliner SN 1991T (Amis(B) « 0.94mag; 
Hamuy et al.|[r996 K 

The host-galaxy reddening, E{B — U)host = 0.206 ± 
0.005 mag measured by SNooPy, is similar to that measured 
by Zheng et al. | 2013[ E{B — V’)host = 0.15 mag) using the 
equivalent-width (EW) of the NalD absorption as a proxy 
for reddening. We adopt our value and the |Cardelli, Clay-| 
ton, fc Mathis] (1989| CCM) reddening law throughout the 


analyses in this work. Our best-fit reddening corresponds to 
a host-galaxy visual extinction of Av.host = 0.64 ± 0.02 mag. 
For the Milky Way reddening, we adopt E{B — U)mw = 
0.14 mag | |Schlegel et al.|[l998| |Schlafly &: Finkbeiner||201l] ) 
and Rv = 3.1. 

SNooPy also determines a distance modulus DM = 
31.49 mag. This gives a luminosity distance of 19.84 Mpc 
for SN 2013dy, which is much greater than the Tully-Fisher 
distance to NGC 7250 of 13.7± 3.0Mpc ( Tully et al.|2009 |. 
However, the distance in SNooPy is derived by fitting the SN 
light curves, so any parameters that affect the light-curve 
properties (e.g., metallicity) could also affect the distance. 
Given that the purpose of this work is to investigate how 
these parameters could affect the SN properties, it would 
be important that we use a SN-independent distance in the 
analysis. Accordingly, we adopt the Tully-Fisher distance to 
NGC 7250 throughout this work. 



0 50 100 150 

Rest-frame days relative to B maximum 

Figure 4. The color curves of SN 2013dy. From top to bottom: 
The B — V, V — R, and V — I color as as the function of phase. 
The filled circles are the data in this work. The color curves from 
SN 2011fe (red open squares) and the [Hsiao et al.| ( [2007^ tem¬ 
plates (blue open triangles) are compared. 


3.2 Color Curves 

The dereddened color curves of SN 2013dy are shown in 
Figures 1^ and For comparison, we also display the color 
curves of SN 2011fe and the Hsiao et al. (20071 templates. 
The Hsiao et al.| | |2007 1 color curves are shifted to match 
SN 2013dy at maximum brightness. For SN 2011fe, we 
adopted the BVRI and JHK photometry studied by [Rich-] 
mond & Smith (20121 and Matheson et al. (20121, respec¬ 


tively. Here the dereddened color curves of SN 2011fe are 
directly compared without applying any shifts. 

We find that after dereddening by our adopted red¬ 
dening parameters, the SN 2013dy and SN 2011fe B — V 
and V — R color curves are generally consistent. However, 
SN 2013dy is slightly bluer (weighted-mean A(B — U) = 
—0.07 ± 0.01 mag and A{V — B) = —0.06 ± 0.01 mag) than 
SN 2011fe at early phases {t ^ 30d). The color curves of 
SN 2013dy also peak at later times (At « 3 days for B — V) 
and have shallower slopes at later phases {t > 30 d) than 
SN 2011fe. Burns et al. (2014| showed that the time of B — V 


maximum is strongly correlated with Ami 5 (B), with faster- 
declining SNe la having earlier B — V maxima and steeper 
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Figure 5. The same as Figure but using V — J and V — H 
color instead. 


late-time B — V slopes than SNe la with slower decline rates. 
Given that SN 2013dy is a slower decliner (Ami 5 (i?) = 
0.9mag) than SN 2011fe (Ami 5 (B) ~ l.lmag), our results 
are consistent with the predicted trends. 

The color curves of SNe 201 Ife and 2013dy have larger 
differences at longer wavelengths. We find that SN 2013dy is 
bluer than SN 201 Ife mV — I a.t all epochs. The same trends 
are also found for the V — J and V — H colors, with even 
greater differences. Finally, relative to SN 2011fe, SN 2013dy 
has much less premaximum V — J and V — H color evolution. 

The color curves of the Hsiao et al. ( 2007| templates are 
generally similar to SN 2013dy at early times. However, they 
fail to reproduce the trends at t > 30 d, where the templates 
are bluer than SN 2013dy. 


3.3 Bolometric Light Curve 

With our exquisite optical/NIR photometry and multi¬ 
epoch UV spectroscopy, we are able to construct a pseudo- 
bolometric light curve which covers 1600-18,000 A, which 
covers essentially the entire SN emission. Since the amount 
of light beyond these wavelengths is likely a very small frac¬ 
tion of all emission, we will refer to the pseudo-bolometric 
light curve as a “bolometric” light curve, but recognize that 
there may be small differences with the true bolometric light 
curve. For all measurements, we correct for Milky Way and 
host-galaxy reddening, using the parameters derived in § |3.1| 

The KAIT BVRI optical light curves, which cover ~ 
3400-9700 A, cover phases from —17d to +340 d. Because 
of the high cadence of the KAIT observations, we do not 
use the RATIR ri light curves to construct the bolometric 
light curve. Integrating the total flux through all filters, we 
construct a BVRI optical light curve. While the bulk of the 
SN emission is in the optical, a significant amount of light 
is emitted in the UV and NIR. 

Using the 10 epochs of HST spectroscopy (covering 
phases of —6.6 to +20.8d), we can construct a UV light 
curve covering ~ 1600-3400 A. The UV flux is calculated 
by integrating the SN spectrum from 1600 A to the effec¬ 


tive wavelength of the B band at that particular epoch. 
For the 10 epochs with UV spectra, we determine the to¬ 
tal UV+optical flux and the corresponding fraction of light 
emitted in the UV. Linearly interpolating the fraction of UV 
light, we estimate the fraction of light emitted in the UV for 
all epochs of the BVRI optical light curve between the first 
and last epochs of UV spectroscopy. 

A NIR light curve was constructed using the YJH 
RATIR light curves. These data cover ~ 9700-18,000 A. 
We did not include the RATIR .Z-band photometry since 
a large portion of the filter overlaps with the KAIT I band. 
The NIR light curve was generated using a procedure simi¬ 
lar to that of the optical light curve. Since our dataset does 
not cover the flux in the K band, we assume a linear de¬ 
cline from the H-band effective wavelength to zero flux at 
18,000 A. Similar to what was done for the UV light curve, 
the fraction of the NIR flux in the optical+NIR light curve 
is calculated. Again, linearly interpolating the NIR fraction, 
the corresponding NIR flux at each individual epoch of the 
optical light curve (for those in the phase range covered by 
the NIR photometry) is determined. 

The final UVOIR bolometric light curve was generated 
by combining the optical light curve with the fraction of 
light determined to be emitted at UV and NIR wavelengths. 
The result is shown in the left panel of Figure Note that 
the bolometric light curve in this work is only valid from 
t ~ — lOd to +22 d, the phases where both HST UV spec¬ 
troscopy and RATIR NIR photometry are available. 

Figure also presents the fraction of the total light 
emitted in the UV and NIR. We find that the UV contribu¬ 
tion peaks at the earliest epochs, with 19% of the bolometric 
light being emitted at UV wavelengths at t « —10 d. The 
fraction of light emitted in the UV declines steeply, becom¬ 
ing only 4% of the total at +21 d. 

The NIR fraction shows a different evolution from the 
UV fraction, being 6% of all light at — lOd, dropping to a 
minimum of 4% at +10 d, and then increasing to 10% at 
+21 d. The NIR contribution overtakes that of the UV at 
~+15 d. 


For comparison, we also show the fractions calculated 
for SN 2011fe (Pereira et al. 20131. The overall trend of 
the UV fraction decreasing and the NIR fraction increasing 
with time also occurs for SN 2011fe. However, SN 2013dy 
emits more light in the UV and less light in the NIR than 
SN 2011fe at all phases. At its peak, SN 2011fe emitted 13% 
of its light in the UV (compared to 19% for SN 2013dy). 
At +20 d, the UV fraction for both SNe is similar (2% and 
4% for SNe 201 Ife and 2013dy, respectively). This trend of 
more UV emission for SN 2013dy with diminishing differ¬ 
ence with time is also seen directly in the UV spectra (see 
§ |4.3[ ). For SN 2011fe, the NIR contributes 9% of the total 
light at —10 d, decreasing some to be 5% at +10 d and then 
increasing to 10% at +20 d, again making SNe 2011fe and 
2013dy be similar at the later phases. 

The bolometric luminosity of SN 2013dy is shown in 
the right panel of Figure |6| Here, we use the Tully-Fisher 
distance to NGC 7250 (13.7 ± 3.0Mpc; |Tully eTaTjl^OOg i. 
Fitting a B-spline to the bolometric light curve, we de¬ 
rive a peak bolometric luminosity of Lboi.peak = lO.Ota'g ^ 
10"^^ erg s“^, where the uncertainty in the peak bolometric 
luminosity is dominated by the uncertainty in the distance 
to SN 2013dy. If the distance derived in SNooPy (19.8 Mpc; 
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Rest-frame days relative to B maximum Rest-frame days relative to B maximum 

Figure 6. Left: The UVOIR bolometric flux of SN 2013dy (open circles). The UV flux fractions determined by HST UV spectroscopy 
are represented by blue filled squares. The NIR flux fractions determined by RATIR NIR spectroscopy are represented by red filled 
triangles. The blue and red dashed line represents the UV and NIR flux fractions of SN 201 Ife | |Pereira et al.|2013[ |, respectively. Right: 
The same as left panel, but using bolometric luminosity instead. The UVOIR bolometric luminosity of SN 2013dy is shown in blue filled 
circles. Here the luminosity distance determined using Tully-Fisher relation (13.7ifc 3.0 Mpc; |Tully et al.|200^ is adopted. The blue curve 
represents the B-spline fit to the data. The grey area shows the Icr range of the bolometric luminosities considering the uncertainties 
from the distance. The bolometric luminosity of SN 201 Ife from[Pereira et al.|||201^ is overplotted (red curve). 


Table 8. The host-galaxy photometry used in this work. 


Table 9. Summary of host-galaxy parameters derived by Z-PEG. 


SDSS u 

SDSS g 

SDSS r 

SDSS i 

SDSS 2 



Lower limit 

Best 

Upper limit 

(mag) 

(mag) 

(mag) 

(mag) 

(mag) 

log Mgteiiar 

(M0) 

8.48 

8.67 

8.84 

13.225(4) 

12.654(2) 

12.352(2) 

12.635(2) 

12.390(4) 

log SFR 

(Moyr-i) 

-0.47 

-0.26 

-0.06 






logsSFR 

(yr 

-8.94 

-8.92 

-8.90 


see § 3.11 is adopted, we derive a peak bolometric luminosity 
of I/boi,peak = 2.09 X 10"*® erg s“^, which is ~ 2 times larger 
than that derived from the Tully-Fisher distance. 

We calculate Amis from the bolometric light curve 
by fitting the light curve with a B-spline. We find that 
SN 2013dy has bolometric Amis = 0.98 mag, which is larger 
than its B-band decline rate Amis(B). Using the same 
method, we determine Amis = 0.99 mag for SN 2011fe. 
This suggests that SN 2013dy and SN 201 Ife have nearly 
identical bolometric decline rates. Zheng et al. I poT^ deter¬ 
mined a rise time (here the rise time is defined as the time 
elapsed from first-light to B-band maximum brightness) of 
~ 17.7 days for SN 2013dy. This value is consistent to that 
of SN 2011fe (17.6 days; Firth et al.|2015 1, and close to the 
average of normal SNe la (17.4 days; |Hayden et al.|2010[ ). We 
find there is hint that SN 2013dy is likely to have longer rise 
time than SN 201 Ife based on the bolometric light-curves 
from t ~ — lOd to maximum brightness. However, the bolo¬ 
metric rise time of SN 2013dy is difficult to constrain directly 
from the bolometric light-curve due to the lack of UV and 
NIR observations in early times. Using the time of first light 


derived by Zheng et al. (20131, we find a bolometric rise time 


for SN 2013dy of ~ 16.47 days. Assuming the Tully-Fisher 
distance of NGC 7250 is correct, this will result in a ®®Ni 
mass estimate of ~ 0.44 Mq for SN 2013dy, using Eq. (6) of 


Stritzinger & Leibundgut (20051. 


3.4 Host Galaxy 

The host galaxy of SN 2013dy, NGC 7250, is a blue, late- 
type (likely between Sc and Im) galaxy. Using SDSS ugriz 


photometry (Ahn et al. 20141 of NGC 7250 (listed in Ta¬ 
ble [8|l, we derive physical parameters for the host galaxy. 
The details of the procedure are described by [Pan et al.| 
( |2014[ ). Here we briefly summarize the process. 

The host-galaxy stellar mass (Mateiiar) and SFR are 
determined using the photometric redshift code Z-PEG (Le 
Borgne fc Rocca-Volmerange|2002 1. Z-PEG fits the observed 
color of the galaxies with galaxy SEDs from 9 different spec¬ 
tral types (SB, Im, Sd, Sc, Sbc, Sb, Sa, SO, and E). We cor¬ 
rect for Milky Way extinction, and allow for a further fore¬ 
ground dust screen varying from B(B — U) = 0 to 0.2 mag 


in steps of 0.02 mag. A Salpeter (19551 initial-mass function 
(IMF) is assumed. 

The best-fit Mgteiiar and SFR determined by Z-PEG are 
log(M/MQ) = 8.67 and 0.55 Mq yr“^, respectively. This in¬ 
dicates that NGC 7250 has a relatively low Msteiiar (e.g.. 


the Milky Way has a log(M/MQ) = 10.79; Flynn et al. 
20061. Based on the mass-metallicity relation derived by 


Tremonti et al. ( 2004[ ), we inferred the gas-phase metallicity 
of NGC 7250 to be 12 + log (0/H) = 8.49. This value is well 
below the solar value (8.69; [Allende Prieto et al.|2001 l, but 
only slightly larger than the metallicity of the Large Mag¬ 
ellanic Cloud (8.40; [Garnett|[l999[ ). Spectroscopic observa¬ 
tions are necessary to precisely constrain the metallicity of 
the host galaxy. Using the Mateiiar derived above, we can 
then determine the SFR relative to the underlying galaxy 
stellar mass, which is called specific SFR (sSFR), the SFR 
per unit Mateiiar ( Guzman et al.|[l997 K Expressing sSFR in 
units of yr”'^, this gives a logsSFR = —8.92 for NGC 7250, 
which is a particularly large value (the Milky Way has a 
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Figure 7. Spectra at t = — 16d and +0d demonstrating the key 
spectral features measured in this work. 


logsSFR = —10.18; Diehl et al.|2006 1 and provides evidence 
that there are strong star-forming regions in the host galaxy. 
A summary of the Z-PEG results is given in Table 


4 SPECTROSCOPIC ANALYSIS 

In this section, we analyze the spectroscopic data. We mea¬ 
sure pseudo-equivalent widths (pEWs) and velocities of key 
spectral features in the UV, optical, and NIR, and study 
their temporal evolution. The variations of narrow absorp¬ 
tion features are also examined. We investigate the spectral 
properties of SN 2013dy by comparing them with those of a 
large sample of SNe la. 


4.1 Photospheric Optical Spectroscopy 


The well-observed spectral sequence of SN 2013dy makes 
it an ideal target to study individual spectral features and 
their temporal evolution. Using the method outlined by 


Maguire et al. (20141 and Pan et al. (20151, we measured 


several spectral features. The key features of interest are 
the pEWs and velocities of the O i triplet, C ii A6580 and 
CiiA7234, SiiiA6355, SiiiA5972, SiiiA4130, and the Can 
NIR triplet. Here we measured Call NIR instead of Call 
H&K since it provides us cleaner measurements of Call ve¬ 
locity and pEW without contamination from other features 


( |Foley| 2013[ [Childress et al.||2014[ [Maguire et al.|[2014| ). In 


Figure 7[we show examples of our photospheric spectra with 

the relevant spectral features marked. Here we are able to de¬ 
compose the high-velocity feature (HVF) and photospheric 
velocity feature (PVF) in the Call NIR triplet (see Figure]^ 
for examples). 

Figure!^ shows the velocities and pEWs of the spectral 
features as a function of phase. Zheng et al. (20131 discov¬ 


ered the presence of both C ii A6580 and C ii A7234 absorp¬ 
tion in the earliest spectra, corresponding to 16 days before 
B maximum brightness, which shows evidence of unburned 
material in the outer layers of the ejecta. In this work we con- 
hrm the detection of C ll and measure physical values similar 


to that of Zheng et al. (20131. The C ii A6580 and C ii A7234 


Rest wavelength (A) 



Rest wavelength (A) 



Rest wavelength (A) 


Figure 8. From top to bottom: The Can NIR absorption at —16, 
+5 and +15 days, respectively. The black curve shows the range 
we fit the line profile. The blue continuum regions near the left 
and right side of the line profile are the regions we select to fit 
the pseudo-continuum. The red solid line is the best fit to the 
line profile, which is the superposition of all the line components 
used for the fitting (represented in dashed lines). The green curve 
represents the superposition of all the components in Can NIR 
photospheric-velocity feature (PVF). The purple curve represents 
the superposition of all the components in the Can NIR high- 
velocity feature (HVF). 
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Figure 9. Top: The SN ejecta velocity as a function of phase for 
various spectral features. Bottom: Same as top panel, but using 
pEWs instead. 


features were found with velocities of 16,300 km s“^ and 
14,800 kms“^ at t = —16d, respectively. However, both fea¬ 
tures quickly faded away after t — —13 d (~3 days after ex¬ 
plosion). The Ol triplet can also be clearly identified in the 
spectral sequence, with an initial velocity of 15,300kms“^ 
at t = — 16d and decreasing to 8,900 km s“^ at maximum 
light. 

The velocity of SillA6355 is ~18,800kms“^ at —16d, 
consistent with the value reported by Zheng et al. (20131. 
The silicon velocity then decreased to 10,200 km s“^ at maxi¬ 
mum light, at which point the velocity became relatively sta¬ 
ble, decreasing to only 9,700 kms“^ at -1-20.8 d. The pEWs 
of Si II A6355 increased slightly from our first epoch until 
t « 5d and then faded away at later times. The SillA4130 
and Si ll A5972 features generally follow a similar evolution 
to that of SillA6355, but with lower velocities and weaker 
absorption. 

The Ca ll NIR triplet can be clearly identified in most of 
our spectra, even extending to very late times. The velocity 
of the Ca ll NIR PAfP and Si ll A6355 have similar evolution, 
but the Call NIR HVF has much larger velocities than the 
Can NIR PVF. For the Can NIR HVF and Can NIR PVF, 
we measure a velocity of 25,130kms“^ and 15,955kms“^, 
respectively, at —16 d. However, the velocity of the Ca ll NIR 
HVF dramatically decreases to ~20,000kms“^ after maxi¬ 
mum brightness. The velocity evolution of Call NIR line is 
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Figure 10. The evolution of the Call NIR triplet with phase. 


consistent to that of normal-velocity (NV; SillA6355 veloc¬ 
ity < 11,800 km s“^ at maximum brightness) SNe la studied 
in Silverman et al. (20151. It is also evident from the pEW 


measurements that the Call NIR HVF is stronger than the 
Call NIR PVF at early times. We find that the ratio of 
pEW(HVF) to pEW(PVF) is ~5.3 at -16 d, quickly de¬ 
creases to unity at t ~ —5 d, and is well below unity at later 
times. 

We also found that the Can NIR absorption showed 
increasingly stronger absorption after maximum brightness. 
We display the evolution of the Call NIR triplet in Fig¬ 
ure which clearly shows the strengthening NIR absorp¬ 
tion with phase. The absorption is still clearly detected at 
t = 131 d. 


4.2 High-Resolution Spectroscopy 


Studies of the circumstellar environment of SNe la are criti¬ 
cal in constraining the progenitor models. By using a series 
of high-resolution optical spectra, recent studies have re¬ 
vealed that at least some SNe la have time-varying narrow 
absorption lines (e.g., Patat et al.||2007 Simon et al.||20d9 


Dilday et al. 20121. Some explanations suggest that these 


time-varying absorptions could be caused by the interaction 
between SN ejecta and the circumstellar medium (CSM), 
which likely originates as gas outflows from the progenitor 
system. 

In this work, we examined the high-resolution spectra 
at wavelengths corresponding to absorption from Ca H&K, 
Na D AA5890, 5896, and K I AA7665, 7699 at roughly the 
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Figure 11. Left: The Na D absorption in the spectrum taken with VLT/UVES. A bar showing the velocity scale is presented on top of 
the spectrum. Right: The EWs as a function of phase of the blue and red components for the Na D1 line as well as the K I A7665 line. 
The spectra taken with different instruments are shown in different colors: HET/HRS (blue), Mercator/HERMES (red), and the rest of 
the instruments listed in Table(black). 




recession velocity of the host galaxy. We also examined ab¬ 
sorption from the 5780 A diffnse interstellar band (DIB). 
The recession velocity of NGC 7250 is sufficiently large 
{cz = 1166kms“^) that there should be no confusion be¬ 
tween interstellar and/or circumstellar gas in NGC 7250 and 
gas in the Milky Way. 


Figure (left) displays the Na D absorption in the 
-1-1.6 day spectrum of SN 2013dy. There are two major ab¬ 
sorption features at roughly —85 and -1-10 km s“^, with the 
bluer component being stronger and saturated. From the 
K I A7665 line, which is not saturated, we determine that 
the —85kms“^ component is composed of two separate ve¬ 
locity components at —83.7 and —88.4kms“^, with the red¬ 
der component being stronger. Using the [Sternberg et ah] 
(20111 classification scheme, where the profile is classified by 
the velocity of absorption features relative to the strongest 
Na absorption, SN 2013dy has a “Redshifted” profile; how¬ 
ever, the bluest (—88.4kms“^) component as seen in K 1 
is weaker than the strongest absorption (—83.7kms“^), 
making the classification somewhat complicated. Using the 


Maguire et al. (20131 classification scheme, where the pro¬ 


file is classified by the velocity of all features relative to the 
systemic velocity, SN 2013dy has a “Blueshifted” profile. 


Since the resolutions for the various instruments differ 
significantly, we focus primarily on the equivalent widths 
(EWs) of the absorption features. However, in the cases 
where we have multiple epochs from the same instrument, we 
also examined the line profiles in detail. Evaluating each sub¬ 
set from the same instrument also reduces cross-instrument 
systematics. There is no obvious change for any line in our 
data. 


We measure the EW of the blue and red components 
of the D1 and D2 lines separately. We do not fit the fea¬ 
tures, but rather simply sum the data over a wavelength 


range covering the features in all spectra. This avoids the 
complication that the higher-resolution spectra have clear 
flat-bottomed profiles, while lower-resolution spectra do not. 
This fiat-bottomed feature has two velocity components (as 
determined by the K I feature) offset by < 5kms“^, and at 
lower resolution they may not be completely resolved. We 
also measure the EW of the entire profile (for D1 and D2 
separately) over a wavelength range covering the entire fea¬ 
ture. These measurements should be sensitive to any changes 
in absorption at velocities different from that of the strong 
absorption features. 


The statistical uncertainties were determined from the 


root-mean-square scatter near the Na lines (e.g., Leonard 
fc Filippenko|200T 1. There are additional systematic uncer¬ 


tainties related to the data reduction (and particularly the 
telluric correction) as well as the continuum determination. 
We measure the latter uncertainty by varying the order and 
spacing of breakpoints when fitting the B-splines to the con¬ 
tinuum, measuring the EW, and setting the standard devi¬ 
ation of the measurements to be the uncertainty. We de¬ 
termine the former uncertainty by comparing our measure¬ 
ments taken with different instruments on the same night 
(with the assumption that lines should not vary significantly 
in the span of a few hours). We present our measurements 
of the EWs for the blue and red components for the Na D1 
line as well as the K I A7665 line in Figure [TT] (right). 


There is no evidence for variable absorption from any 
of the narrow absorption features in our data. Although the 
scatter for a given measurement is significantly larger than 
the uncertainties in the measurements (with ~ 10), there 
is no coherent change with phase. Instead, we believe that 
a systematic bias affects the HET/HRS observations, which 
have systematically high EWs for all components of Na D. 
While we do not attempt to correct this potential bias in the 
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Figure 12. Left: The comparison of SN 2013dy (black) and SN 2011fe (red) spectra under 5000 A. The maximum-light spectra of 
SN 2011by (green) and SN 2011iv (blue) are also compared. All spectra are dereddened and normalized accordingly (the normalized 
region is 3500-5000 A). Right: The flux-ratio spectra (llfe/13dy) are shown in black. The predicted flux-ratio spectra with differences 
in metallicity factor f are overplotted. Here the models from Days 7, 10, and 15 of [Lentz et al.| ( |2000^ are compared to the dataset at 
t = —7, —3, and 0 days relative to the R-band maximum brightness, respectively. The grey dashed line represents the line of equality. 


measurements, we note that the scatter is typically < 10% of 
the measurement, indicating that the potential bias is likely 
small. The peak-to-peak variations up to ~ 40 mA are also 
expected owing to the fractal patchy structure of the ISM 
| |Patat et al.|20To| . 

In our highest S/N spectra, we detect the 5780 A DIB 
feature, which is consistent with a single component having 
a velocity of —20± 10kms“^. This component is marginally 
consistent with being at the same velocity as the redder 
-|-10kms“^ Na D absorption co mponent. The DIB h as an 
EW of 60 ± 10 mA. Using the Phillips et al. (20131 rela¬ 
tion, this corresponds to Ay = 0.3 ± 0.2 mag. This esti¬ 
mate is comparable to, though lower than, the extinction 
measurement determined from the SN 2013dy photometry 
(Av = 0.64 mag; see § 3.11. 


4.3 UV Spectroscopy 


We obtained 10 epochs of HST/STIS spectra in this work. 
The multi-epoch UV spectra enable an investigation of the 
UV spectral evolution at wavelengths as short as 1600 A. In 
Figure [T2| we compare the HST spectra of SN 2013dy with 
those of SN 201 Ife at similar epochs. We also compare to 
SNe 2011by and 2011iv at maximum brightness. SNe 2011by 
and 2011fe have nearly identical optical colors, light-curve 


shapes, and spectra, but have different UV continua (Fo- 


ley fc Kirshner|[2013| Graham et al.|[2015a |. SN 2011iv is a 


spectroscopically normal SN la, but presents a relatively fast 
decline rate (Ami 5 (B) « 1.7mag; Foley et al. 2012bI. All 
spectra were dereddened and normalized to have the same 
mean flux between 3500 and 5000 A. 


At all epochs except for t = +21 d, SN 2013dy has 
excess flux at 2800 < A < 4000 A relative to SN 2011fe. 
Examining this wavelength range at maximum brightness, 
we find a trend between the amount of flux in this region 
(relative to the rest of the spectrum) and light-curve shape. 
Specifically, the continuum in this region is strongest for 
SN 2013dy {Amis{B) = 0.92 mag), followed by SNe 2011by 
and 2011fe {Ami^{B) ~ 1.1 mag), and finally SN 2011iv 
(Ami 5 (B) « 1.7mag). Therefore, the amount of flux in this 
region may be linked to the amount of ®®Ni generated in 
the explosion. We also note that the maximum-light spectra 
of SNe 201 Iby and 2013dy are relatively similar for A < 
2300 A, with both having relatively less flux than SN 201 Ife 
at these wavelengths. SN 2011fe also has excess flux at these 
wavelengths, relative to SN 2013dy, at earlier epochs. 


We further examine the UV evolution of SN 2013dy 
using the flux-ratio spectra (fisdy/fiife) shown in the right- 
hand panel of Figure Followed the method in [Foley fc] 


Kirshner (20131, we overplot the flnx-ratio spectra predicted 
by the model of |Lentz et al.| | |2000[ ) with differences in metal¬ 
licity factor We find the SN 2011fe/SN 2013dy flux-ratio 
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Figure 13. Same as the right panel of Figure but using the 
flux-ratio spectrum of SN 2011by and SN 2013dy instead. 


spectra are better matched by the Lentz et al. (20001 flux- 
ratio spectra with >30 times metallicity difference. This 
might suggest that SN 2013dy has a higher progenitor metal¬ 
licity than SN 2011fe; however, the other differences between 
these SNe (light-curve shape and spectra; see § 4.5 for de¬ 
tails) complicate this possibility. 

For comparison, we also determine the 
SN 2011by/SN 2013dy flux-ratio spectrum at maximum 
light (see Figure 131. The result shows that SN 2013dy 
and SN 2011by have similar continua in the UV. This 
is consistent with the conclusion from |Foley fc KirshriSI 
(20131 and Graham et al. (2015a I that SN 2011by appears 
to have a higher progenitor metallicity than SN 2011fe. 

In addition to the difference in the continuum, the po¬ 
sitions of the UV spectral features of SN 2013dy are system¬ 
atically bluer than those in th e spectra of SN 2011fe . While 
lacking data at A < 2900 A, Maguire et al. (20121 saw a 
similar trend for the features at 2900 < A < 3300 A when 
comparing the mean spectra of SNe la with slower decline 
rates to those having faster decline rates. 


Foley et al. (20081 defined the “UV ratio” as i?uv = 
f\ (2770 A)//a (2900 A), finding that the UV ratio is strongly 
correlated with the luminosities of SNe la. Brighter SNe la 
tend to have smaller Ruv than their fainter counterparts. 
For SNe 2013dy and 2011fe, we measure i?uv = 0.36 and 
0.42, consistent with their light-curve shapes. 

The strong absorption feature at ~3000 A (which is at¬ 
tributed to Fell A3250; Branch fc Venkatakrishna|1986| ) can 
be clearly identified in all UV spectra. However, the blending 
with many IGE lines in Fe ll A3250 complicates the velocity 
measurement. Here we only measure the pEW of the feature 
as a whole at each epoch and present the result in Figure]^ 
The pEW of Fe ll A3250 increases slightly from ~ 10 A at 
— 7d to ~60A at -1-21 d. Our results compare well with the 
evolution studied by Foley et al. (20081 using a larger sample 
of SNe la. 


4.4 Nebular Spectra 

By ~100d after maximum brightness, SNe la enter the so- 
called ‘nebular phase.’ At this point, the SN becomes opti¬ 
cally thin and photons can escape from the very center of 
the ejecta. Thus, observations at these epochs are strongly 
constraining for SN la explosion models. For these phases, 
the SN spectrum is dominated by forbidden emission lines of 
IGEs (e.g., iron, nickel and cobalt). In this work, we obtained 



Rest-frame days relative to B maximum 


Figure 14. The velocity determined from the emission lines 
in SN 2013dy nebular spectra as a function of phase. The ve¬ 
locity below/above the dashed line means the emission line is 
blueshifted/redshifted relative to the rest wavelength. 
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Figure 15. Late-time spectra of SN 2013dy (black) and 
SN 201 Ife (red). The spectra are dereddened and normalized ac¬ 
cordingly (the normalized region is 5000-6000 A). 


6 nebular spectra at f > 100 d (three with 100 < t < 150 d 
with the remaining at f > 300 d). 

For each spectrum, we measured the velocity shift of 
the [Feiii] A4701, [Fell] A7155, [Goiii] A5891 and [Niii] A7378 
emission line features in the nebular spectrum. Figure 
shows the velocity of emission lines as a function of phase. 
[Feiii] A4701 is one of the strongest spectral features in the 
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nebular phase, and we are able to measure its velocity in 
most of the nebular spectra. It has a blushifted velocity of ~ 
— 2600 km s“^ relative to the rest wavelength of [Fe ill] A4701 
at t ~ food. The velocity then decreases to —720 kms“^ at 
t = 333d, and is redshifted to ~1300 kms“^ at t = 423d. 


§ 3.31 may be underestimated. Given the relative light-curve 


Similar trend was also found in Silverman et al. (20131 with 
a larger sample of SNe la. 

[Com] A5891 shows a similar evolution to [Feiii] A4701, 
but has a velocity that is ~ 1000 kms“^ closer to zero 
shift. [Fell]A7155 and [Nill]A7378 are blended in wave¬ 
length space and can only be clearly decomposed in the 
spectra at t = 333 d and 423 d. We deblend the line fea¬ 
ture by fitting the line with a double Gaussian prohle. The 
velocities of [Fell] A7155 and [Nill] A7378 are similar to that 
of [Felll]A4701 at t = 333 d, although an opposite evolu¬ 
tion is found afterwards (the velocities of [Fell] A7155 and 
[Nill] A7378 are further blueshifted at 423 d). 

In Figure we compare the late-time spectra of 
SN 2013dy (starting ~1 month after maximum brightness) 
with similar phase spectra of SN 201 Ife. The results show 
that both SNe are very similar through t « 100 d. How¬ 
ever, there are noticeable deviations for the [Fell]A7155 
and [Nill]A7378 features as well as the spectral region 
around 8500-9000 A (dominated by [Fe ll] lines) at t = 333 d. 
SN 2013dy appears to have a relatively weaker [Fe ll] A7155-I- 
[Nill]A7378 feature than that of SN 2011fe. At this late 
phase, the [Niii[ A7378 pEW of SN 2011fe (1004 ± 2 A) is 
about 1.6 times larger that of SN 2013dy (630 ± 6 A), but 
both SNe have similar [Fell] A7155 pEWs (739 ± SA and 
892 ± 1 A for SN 2013dy and SN 2011fe, respectively). 


Following the diagnostic in Graham et al. (2015aI, 


we determine the relative strength of [Ni ll] A7378 to 
[Fell] A7155, but using pEW instead of the flux of the emis¬ 
sion line. Here the pseudo-continuum used to calculate the 
pEW of emission line is dehned with the local minima on 
each side of the line profile. We measure a pEW ratio of 0.85 
and 1.12 for SN 2013dy and SN 2011fe, respectively. Since 
the radioactive ®®Ni has a half-life of ~6d, all of the nickel 
emission in the nebular spectra at t > 300 d is produced by 
stable ®®Ni. However, the [Fell] A7155 emission is produced 
by the combination of stable iron and radioactive ®®Ni. This 
degeneracy complicates our analysis. 

For two identical SNe la, except for the amount of 
(stable) ®‘^Fe, radioactive ®®Ni, and (stable) ®®Ni gener¬ 
ated in the explosion, their nebular spectra would be sim¬ 
ilar, with the largest discrepancy expected to be the flux 
ratio of [Niii] A7378/[Feii] A7155. If the amount of stable 
iron-group elements is the same, but more ®®Ni is gener¬ 
ated, which will mostly decay to ®®Fe at the time a neb¬ 
ular spectrum is observed, one would naively expect the 
[Ni ii] A7378/[Fe ii] A7155 flux ratio to be smaller than in 
an object with smaller ®®Ni. Similarly, if the amount of 
®®Ni is the same, but the amount of stable iron-group 
elements is smaller, then one would expect a smaller 
[Nill] A7378/[Feii] A7155 flux ratio. 

Since SN 2013dy has a smaller [Nill] A7378/[Feii] A7155 
flux ratio than that of SN 2011fe, SN 2013dy naively had a 
smaller amount of stable iron-group elements or a larger 
amount of ®®Ni (or some combination) than did SN 201 Ife. 
For the latter case, and given the correlation between SN la 


shapes for the two SNe, with SN 2013dy having a broader 
light curve than SN 2011fe, this is a reasonable explanation. 
A separate distance measurement to the host galaxy (e.g., 
using Gepheid distance) is necessary to verify our results. 

For the former case (SN 2013dy having a smaller 
amount of stable iron-group elements than SN 201 Ife), the 
relative [Nill] A7378/[Fell] A7155 flux ratio would suggest 
that SN 2013dy produced a smaller amount of stable to 
radioactive iron-group elements, consistent with SN 2013dy 
having a smaZ/er metallicity than SN 2011fe, contrasting the 


results from our UV data (see Section 4.31. 

Again, the real situation may be more complicated, and 
some degeneracies need to be considered here. SNe la with 
higher metallicities not only produce more ®®Ni, but also 
more stable iron ( [Timmes et al.|[2003 l. That means the ra¬ 
tio [Nill] A7378/[Fell] A7155 may not be necessarily larger 
with higher metallicity. In addition, the prediction made by 
Timmes et al.| ( |2003| ) is for the exact same explosion ener¬ 
gies. As we will discuss in § |4.5[ this is unlikely to be the 
case for SN 2013dy and SN 2011fe. Physical conditions other 
than progenitor metallicity could also contribute here (e.g., 
central density; ILesaffre et al.||2006[ [Krueger et al.|[20T2l 


Seitenzahl et al. 12013 1. Seitenzahl et al. (20131 showed that 


in Chandrasekhar-mass models the influence of central den¬ 
sity (at ignition) on the ratio of [Nill] A7378/[Feii] A7155 is 
even greater than that of metallicity. They found that the 
models with higher central density generally have a higher 
ratio of [Nill] A7378/[Fell] A7155 (assuming they have equal 
metallicity). It is currently difficult to disentangle the two 
situations for singular objects. Nevertheless, the data ap¬ 
pear to be more consistent with SN 2013dy having produced 
more ®®Ni than SN 201 Ife, with metallicity being less of an 
influence. 


4.5 Spectroscopic Subclasses 

Previous studies which examined large sample of SNe la 
and their spectra at maximum light demonstrated the ex¬ 
istence of several subclasses of SNe la (Benetti et al.|[^005| 


Branch et al. 20061. In Figure 16 (left) we show the sub¬ 


classes of SNe la as determined by the gradient of Si ll A6355 
velocity and B-band decline rate ([Benetti et al[] 

20051. The SN la sample in Benetti et al. ( 2005|) is split into 


three groups: a ‘low-velocity gradient’ (LVG) group, which 
have low SillA6355 velocity gradients and natrually have 
lower values of Ami 5 (B); a ‘high-velocity gradient’ (HVG) 
group, which have high SillA6355 velocity gradients and 
low Ami 5 (B); and a ‘FAINT’ group including SN 1991bg- 
like events that have higher values of Ami 5 (B) and nat¬ 
urally have larger velocity gradients. SN 2013dy shows a 
low SillA6355 velocity gradient (—15kms“^ day“^; mea¬ 
sured from maximum-light to -1-21 d) and a slow decline rate 
(Ami 5 (B) = 0.92mag). Using these spectral indicators, we 
can classify SN 2013dy as a LVG SN la. 


Figure 16 (right) shows the subclasses of SNe la ac¬ 


peak bolometric luminosity and the mass of °°Ni (Arnett 


1979 19821, the Tully-Fisher distance to NGC 7250 (see 


cording to the pEW of Si ll A6355 and Si ll A5972 (Branch 
|et al.|[2006[ ). [Branch et al.| ( |2006| ) split the SN la sample 
into four groups: A ‘shallow-silicon’ group, which have low 
pEWs for both SillA6355 and SillA5972 lines (and include 
SN 1991T-like objects); a ‘core-normal’ group, which have 
homogeneous and intermediate pEWs; a ‘cool’ group (simi- 
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Figure 16. Left: The gradient of SillA6355 velocity (kms“^ day“^) as the function of S-band decline rate Ami 5 (_B). The sample 
and subclasses studied in |Benetti et al.1 ( |2005^ is overplotted. The sample is split into three different subclasses: ‘low-velocity gradient’ 
(LVG; open squares), ‘high-velocity gradient’ (HVG; open circles) and ‘FAINT’ (open triangles). Right: SillA5972 pEW as the function 
of Si IIA6355 pEW. The sample and subclasses studied in [Branch et al.| ( [2009^ is overplotted. The sample is split into four different 
subclasses: ‘shallow silicon’ (open triangles), ‘core normal’ (open circles), ‘cool’ (open downward triangles) and ‘broad line’ (open squares). 
The position of SN 2013dy is represented by the red filled star. We also mark the positions of a variety of SNe la for comparison. 


lar to the FAINT group in Benetti et al.|[2005 l with strong 
SiiiA5972 lines relative to SiiiA6355; and a ‘broad-line’ 
group, which present strong SillA6355 absorptions (and 
have much overlap with the HVG group). In this scheme, we 
classify SN 2013dy as a ‘shallow-silicon’ SN la corresponding 
to its relatively weak Si ll A6355 and Si ll A5972 absorptions. 


In addition to the samples of Benetti et al. (20051 
and [Branch et al.j ( |2006[ ), we specifically added SNe 2009ig, 
2011by, 2011fe, 2011iv, and 2014J to F igure [16| All of these 
SNe have maximum-light UV spectra (Foley et al. ]|2012a|b[ 


Foley fc Kirshne'r|[2013[ [Foley[|^13[ [Foley et al.[[2014| [Maz-[ 


zali et al.||2014 K While SNe 2009ig and 2014J have similar 

decline rates as SN 2013dy (Ami5(i3) = 0.89 and 0.95mag, 
respectively), both have larger Sill A6355 velocity gradients, 
with SN 2014J having a velocity gradient somewhat inter¬ 
mediate to the LVG and HVG classes. Nonetheless, all three 
SNe are classified as LVG objects. SN 201 Ife, on the other 
hand, is intermediate to the LVG and HVG groups. 

Using purely spectral parameters, we classify 
SNe 2011by and 2011fe as Core-normal, SN 2011iv as 
Cool, SN 2014J as Broad-lined, and SN 2009ig as inter¬ 
mediate between Shallow-silicon and Core-normal. This 
indicates that there is no perfect match for SN 2013dy 
among the SNe la with maximum-light UV spectra. How¬ 
ever, from Figure [T^ we see that we are consistently filling 
out the parameter space of SNe la with objects that have 
UV spectra. To this end, SN 2013dy represents an extreme 
object in velocity gradient, Ami 5 (B), pEW(SillA6355), 
and pEW(SiiiA5972). 


5 COMPARISON TO MODELS 

With an excellent UV through NIR spectral sequence and a 
well-sampled bolometric light curve, SN 2013dy is an ideal 
object to which one can compare models. Through these 
comparisons, we can constrain the progenitor and explosion 
models of SN 2013dy. For this purpose, we select two pri¬ 
mary models: The one-dimensional deflagration model W7 


I Nomoto et al. 19841 and the three-dimensional delayed- 
detonation N5 model ( [Seitenzahl et aL][2013[ ). These mod¬ 
els are chosen based on their similar bolometric light-curve 
shapes to that of SN 2013dy (e.g., see Sim et al. 2013| ). Spec¬ 
tral sequences were generated for these models using the ra¬ 
diative transfer codes developed by Kromer & Sim (ARTIS; 


20091. The corresponding model spectra were described in 
Krorner fc Sim[ ( f2009[ ) for the W7 model and jSim et al.| ( [?013[ ) 


for the N5 model. In addition to the nominal W7 model, we 
modify the W7 abundance profile to match solar metallicity 

J Asplund et al.[[2009 1 in the outer layers (W7 Zq; see also 
bley et al.[[2012b[). This model is useful for determining if 


the progenitor star had nonzero metallicity. 


In Figure we compare the bolometric light curve of 
SN 2013dy with the model bolometric light curves. Because 
of the large distance uncertainty of SN 2013dy, we shift the 
model light curves such that their flux at the time of B- 
band maximum brightness match that of SN 2013dy. By 
doing this, we are examining the shape of the bolometric 
light curve instead of the absolute luminosity. Future obser¬ 
vations may more accurately determine the peak luminosity 
of SN 2013dy and thus better constrain the explosion mod¬ 
els. 


At early phases (t Od), the model light curves gener¬ 
ated from W7 Zq and N5 models are nearly identical, and 
both are a better match to the data than the W7 model. 
At later phases {t > Od), the N5 light curves start to de¬ 
viate from the data, having a slower decline, while the W7 
and W7 Zq evolve in a similar fashion as SN 2013dy. The 
slower decline found in N5 model could be a reflection of 
the larger iron-group element mass compared to W7 model, 
which will result in a higher opacity and temperature. The 
difference between the W7 and W7 Zq models are negligi¬ 
ble at later epochs since the W7 and W7 Zq models only 
differ in the outer layers. The postmaximum evolution of 
both models should resemble to each other as expected. We 
further demonstrate this effect in Figurej^by investigating 
the details of the spectral sequence. 
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Previous theoretical studies have shown that the UV 
SED is extremely sensitive to the progenitor metallicity as 
it significantly alters the UV opacity, which is dominated 
by line blanketing i Hoeflich et al.|fl998 Lentz et al.||2000 |. 
SNe la with identical optical spectra could have (and have 
shown) large differences in the UV I Foley fc Kirshner|2013 1. 
In Figure |18[ we compare the model spectra with those of 
SN 2013dy. To be consistent, the model spectra are scaled 
with the same multiplicative factor used to shift the model 
bolometric light curve in Figure We exclusively examine 
the HST spectra, since almost all of the discriminating power 
is in the UV. 

For the phases earlier than a week after maximum 
brightness, the W7, W7 Zq, and N5 models have similar 
optical spectra with all generally agreeing with the data. 
However, as already noticed by Foley et al. (2012b i, the W7 
and W7 Zq models start to have large differences at far-UV 
wavelengths (A ^ 2800 A), which is exactly the region shown 


to be sensitive to the progenitor metallicity (e.g., Lentz et al. 
|2000[ ). While the UV flux of the W7 model is significantly 
larger than that of SN 2013dy (>50 times at some wave¬ 
lengths), the W7 Zq and N5 models are generally in agree¬ 
ment with the data. This result seems to be a conflict with 
our previous conclusion since the N5 model does not have so¬ 
lar metallicity admixture like the W7 Zq model. However, in 
the three-dimensional delayed-detonation models, the IGE- 
rich deflagration ashes will rise due to the turbulent motions 
(e.g., Sim et aL]|2013 1. This will lead to a mixture of burn¬ 
ing ashes at higher velocities and cause the suppression of 
UV flux, which is similar to the effect of progenitor metal¬ 
licity in W7 Zq model. For our earliest epoch {t = —6.6 d), 
even the N5 model has more UV flux than the data, but at 
this time, the W7 Zq model is relatively consistent with the 
observations. 

For the phases later than a week after maximum bright¬ 
ness, the models are less consistent with the observations. 
These discrepancies might be caused by systematic effects 
in the model spectra related to approximations in the treat¬ 
ment of excitation and ionisation in our radiative transfer 
simulations. Since full non-LTE calculations are currently 
unavailable for these explosion models, this limitation is ex¬ 
pected to introduce some uncertainty to the radiative trans¬ 
fer simulations in this work. 

To summarize, we find that the W7 Zq model is in good 
agreement with SN 2013dy. In terms of bolometric behavior 
and early-time UV spectra, it performs significantly better 
than the zero-metallicity W7 model. For the earliest epoch of 
spectroscopy and at later times, the W7 Zq model also per¬ 
forms slightly better than the delayed-detonation N5 model. 
The UV data were critical in making this determination. 


6 CONCLUSIONS 

SN 2013dy was discovered only 2.4 hr after explosion in 
NGC 7250, which makes it the SN la with the earliest known 
detection. In this work we present high-quality photometric 
and spectroscopic data. The photometric observations were 
taken from f = —16 to -1-337 d (relative to maximum bright¬ 
ness). Ten epochs of UV through NIR spectra were taken 
with i/ST/STIS. For phases of f = —16 to -1-480 additional 



Rest-frame days relative to B maximum 


Figure 17. The UVOIR bolometric light curve of SN 2013dy 
(same as Figure]^ is shown in open circles. We compare the data 
to model light curves generated from W7 model (blue curve), 
solar-metallicity polluted W7 Zq model (red curve) and the 
three-dimensional delayed-detonation N5 model (green curve), re¬ 
spectively. Here the model light curve shown for N5 model is 
angle-averaged. The model light curves are shifted to match the 
bolometric flux of SN 2013dy at the time of S-band maximum. 


low- and high-resolution optical spectra were obtained with 
a variety of ground-based facilities. 

Our photometric analysis shows that SN 2013dy has 
a relatively slow decline rate (Ami^{B) = 0.92mag) and 
is intrinsically bluer than SN 2011fe. Gonsistent with pre¬ 
vious results, SN 2013dy, being a slow-decling SN la, has 
color curves that peak later and have shallower later-time 
slopes than those of faster-declining SNe la. We gener¬ 
ated a UVOIR bolometric light curve for SN 2013dy (1600- 
18,000 A) with both photometric and spectroscopic datasets, 
but large uncertainties in the distance to SN 2013dy pre¬ 
vented a precise measurement of its peak bolometric lumi¬ 
nosity. 

We also determined the physical parameters of 
NGG 7250, the host galaxy of SN 2013dy. NGG 7250 is 
blue with evidence of strong star-formation. The relatively 
low Msteiiar suggBsts that the environment of SN 2013dy is 
likely metal poor. This seems to be a conflict with the re¬ 
sults suggested by the UV spectroscopy. However, using the 
global properties of host galaxy as proxy to constrain the 
SN progenitor could present large uncertainties. The host 
parameters determined in this work may not well reflect the 
real physical conditions of the SN birthplace. 

In the earliest spectra (1.6 to 5d after the explosion), 
there are strong G II absorption features, indicating the 
evidence of unburned material exists at the outer layers 
of the ejecta. Using the classification schemes defined by 


Benetti et al. (20051 and Branch et al. (20091, we further 


classified SN 2013dy as a SN la with low-velocity gradient 
(LVG) and shallow-silicon absorption. Therefore, SN 2013dy 
is spectroscopically different from Gore-normal SNe similar 
to SN 2011fe. 


SN 2013dy also has strong Call NIR HVF at early 
phases. The velocity of the Call NIR HVF is ~10,000 kms“^ 
higher than that of the Call NIR PVF. We find that the 
Call NIR HVF was stronger than the Can NIR PVF a 
few days after the explosion. However, the Call NIR HVF 
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Figure 18. The HST spectra at different epochs (black curves). We compare the data to model spectra generated from W7 model (blue 
curves), solar-metallicity polluted W7 Zq model (red curves) and the three-dimensional delayed-detonation N5 model (green curves), 
respectively. The model spectra are scaled with the same multiplicative factor used to shift the model bolometric light curve in Figure [TT] 


quickly became weaker than Ca ll NIR PVF immediately af¬ 
ter t = —5d. Although the velocity of the Call NIR triplet 
decreased quickly with phase, we saw increasingly stronger 
absorptions for the Ca ll NIR triplet from maximum bright¬ 
ness until t = -1-131 d. 

We obtained 15 epochs of high-resolution spectra, which 
is one of largest known single set of high-resolution spectra 
for a SN la. Examining narrow absorption features including 
Na I D, K I A7665, and the 5780 A DIB, we find no evidence 
of temporal variability. 

We compared the HST UV spectra to those of 
SNe 2011by, 2011fe, and 2011iv at similar phases. At maxi¬ 
mum brightness, there appears to be a correlation between 
the amount of UV flux at 2800 < A < 4000 A and decline 
rate with SN 2013dy having both the slowest decline rate 
and highest UV flux of the sample. By examining the flux- 
ratio spectra, we find that SN 2013dy has a depressed far-UV 
continuum (A < 2300 A) relative to that of SN 2011fe, con¬ 
sistent with predictions of a higher metallicity progenitor 
for SN 2013dy. However, differences in the two SNe com¬ 
plicate this interpretation. SN 2013dy also shows systemat¬ 
ically blueshifted UV features relative to SN 201 Ife. 

The late-time spectra of SNe 2013dy and 2011fe are 
very similar from t = 40 to 100 d. However, we find that 
SN 2013dy has a smaller ratio of [Niii] A7378/[Feii] A7155 


than SN 2011fe at t ~ 333 d. This may suggest that 
SN 2013dy synthesized more ®®Ni than SN 2011fe assum¬ 
ing both SNe generated the same amount of stable iron and 
®®Ni, but we note that metallicity could have some effect 
here. 


By comparing the bolometric light curve and spectral 
sequence of SN 2013dy to that of SN la explosion models, 
we find that the W7 Zq model (solar metallicity polluted 
W7 model) is in good agreement with the entire SN 2013dy 
dataset. The W7 Zq model is much more consistent than 
the zero-metallicity W7 model and slightly more consistent 
than the delayed-detonation N5 model. SN 2013dy again 
demonstrates that UV data are critical for understanding 
the progenitor metallicity for SNe la. 


SN 2013dy is one of the best-studied SNe la. It is very 
nearby and bright, and was discovered immediately after ex¬ 
plosion. Follow-up observations were obtained consistently, 
spanning both a large time frame and wavelength range. 
This exquisite dataset shows that a detailed study of a single 
well-observed SN la can provide unique information about 
the progenitors and explosions of SNe la. 
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